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Abstract 

In this paper, we show how the non-holonomic control technique can 
be employed to build completely controlled quantum devices. Examples 
of such controlled structures are provided. 

1 Introduction 

In the previous article, we showed that quantum systems become "non-holonomic' 
when perturbed in a certain time-dependent way : as a result of the perturba- 
tion, all global constraints on the dynamics are removed and the system becomes 
fully controlled. The straightforward application of the non-holonomic control 
to groups of interacting two-level systems, the so-called qubits, is a promising 
way to achieve computational algorithms widely discussed in quantum infor- 
matics P 121 El 0- Indeed, quantum computations [S] are achieved through 
combining well-chosen quantum gates [Hj, which are the analog of the classi- 
cal logic gates and consist in particular evolutions of the system chosen as a 
computer. In this context, the non-holonomic control is a precious means for 
implementing any gate and thus performing any computation on an arbitrary 
quantum system. However, the direct control of an A^-qubit computer's evolu- 
tion is a very heavy (or even intractable) computational task when A'^ is a large 
number, since it requires the control of 4^ physical parameters. To overcome 
this impediment one can combine completely controlled cells in a manner which 
depends on the kind of the problem to be solved : the adaptation of the struc- 
ture of the quantum device, obtained as a particular arrangement of controlled 
elements, to the computation to be performed allows one to decrease the number 
of the free control parameters needed, and thus the complexity of the control 
problem. The complete control of the whole compound device is thus assured 
by the controllability of the individual cells as well as their connections with 
each other within the architecture of the computer. 

In the first section of this paper, we provide an example of a unit cell, which 
is completely controlled through non-holonomic interactions. In the second sec- 
tion, we propose two different devices composed of such cells, the arrangements 
of which suit particularly well universal quantum computations and simulation 
of quantum field dynamics, respectively. In the third section, we finally describe 
a toy device that can perform quantum computations on 9 qubits and show in 
particular how it can perform the discrete Fourier transform on 9 qubits. 

2 Completely Controlled Unit Cell 

One way to construct a completely controlled but not immediately universal 
quantum device is to build it up from small parts, called "unit cells", each of 
which is non-holonomic and therefore directly and universally controllable. The 
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proper functioning of the device relies then on the appropriate connection of 
the ceUs. In this way the universahty of the device is obtained indirectly, not 
by applying a huge number of controls, but by smartly connecting the cells and 
choosing the sequence of operations performed. There is no general prescription 
on how to construct a particular device; this requires expertise in the art of 
"programming" the operations of the cells and their interactions. 

2.1 Cell structure 

Fig- n shows an example of a completely controlled unit cell composed of three 
two-level atoms, each with ground and excited states |0) and having distinct 
transition frequencies wj, Wj, and ui^. The atoms are subject to dipole-dipole 
interactions and are coupled to two external fields: an electromagnetic field 
Eu! = SujCOSLot of nearly resonant frequency oj, and a static electric field Eg- 
The dipole-dipole interaction is fixed and determines the principal, unperturbed 
Hamiltonian of the system, Hq , while the external fields provide two controllable 
perturbations, P^^ and Ps- The Hilbert space of the system has a "computational 
basis" of = 2^^ = 8 states, \x) = \x2Xixq) = \x2)\xi)\xo) , x — 0,1,..., 7, 
where the state of the ith atom encodes the ith binary digit oix = X)r= 
a qubit [see Fig.^b)]. The crucial requirement is the non-holonomic character 
of the interaction : H^, P^, Pg, and their commutators of all orders must span 
the linear space of 8 x 8 Hermitian matrices. This is indeed the case for the 
system shown in Fig. ^ for which principal Hamiltonian and perturbations are 
given, in the computational basis and assuming resonant approximation, by the 
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Here Dij = didj / R^^ denotes the dipole-dipolc coupling of the ith and jth atoms 
at distance i?^, with di the ith atom dipole matrix element, and Vi = Sujdi 
is the dipole coupling of the ith atom to the external electromagnetic field. 
The excitation energy dctunings of single atoms Ai = h{ijjf — uj) determine 
the detunings of pairs of atoms Aij = Ai + Aj and the total detuning A<j = 
Ai+ A2 + A3. Their values can be changed by variation of a static electric field 
Es (Stark effect), which results in energy shifts Ai ^ Ai + Aj for single atoms, 
where Aj = aiEs depend on atom-specific electric permeability constants ai, 
and similar shifts Aij = Aj + Aj and A^ = Ai + A2 + A3 for two and three 
atomic dctunings respectively. 

Note that, by a proper choice of Aj and ui, one can set two of three Ai to 
zero. Moreover, for clarity we also set to zero the third Ai, which remains just 
a part of Hq otherwise. Hence, hereafter, all Aj denote just the deviations from 
zero resulting from the variation Es of the Stark field. The last together with 
the amplitude = serve as a time dependent control parameters Cs and Co, 
respectively. The matrices P5 and contain therefore only the permeabilities 
ttj and the dipole moments di respectively. 
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Figure 1: Realization of a unit cell: A compound system of three two-level 
atoms interacting with external electromagnetic and static electric fields, (a) 
The ith atom has ground and excited states |0)i and jl); with excitation energy 
Ai + Ai that can be modified by the static field; transition amplitude in the 
electromagnetic field is Vf, the dipole-dipole coupling of the zth and jth atoms 
is Dij . (b) The computational basis states and their relation to matrix elements 
of the principal Hamiltonian Hq and the perturbations P^^ and Pg of Eqs. Hll3f) . 
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2.2 Cell control 



To exert direct universal control over the unit cell we employ the non-holonomic 
control technique presented in the previous article, (i) We fix N'^ = 64 consec- 
utive time intervals of equal length t = during which the two perturbations 
are alternately applied to the system : in the fcth interval the perturbation is 
Pk = Ps for odd k and Pk — Pui for even k, where k — 1,2, ... , 64. The strength 
of Pk is denoted by Cfc, and corresponds either to or Es, depending on the 
parity of k. Thus, the evolution of the system is governed by a Hamiltonian 
which is constant on each interval: 



H{t) = Ho + CkPk te[{k^ 1)t, kr]. 
(ii) We look for the 64-dimcnsional vector C such that 



(4) 
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To this end, we first solve the "8th root" of Eq. (O, 
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by minimizing the functional Y!fj=o 1% ({cfc}fc=i...8)|^ to 2, where {aj} denote 
the coefficients of the characteristic polynomial of the matrix product in Eq. © . 
This provides a sequence of eight values, ci,C2, . . . ,C8, the repetition of which 
yields the desired vector C . (iii) Finally, we compute the required vector C 
as described in the previous article. If the target evolution UarUtrary — = 
exp(— zTie) is close to the identity (i.e. e is small), one determines the variations 
5Ck to first order in e by solving the linear equations 



64 „p 

fc=i 



(7) 



and, replacing C^^^ by C^"^ + 5C , one repeats the same operation, and so on, 
until one gets the vector C which checks C/ ^ C ^ =Ue with desired accuracy. To 

perform an arbitrary unitary transformation UarUtrary = , with e taking any 
value in [0, 27r] and not necessarily small, we divide the work into "small" steps : 

we apply the transformation U ^c(^)^ = [/<:/„* = {UarUtrary^ " repeatedly 
n* times, with n* determined as described in the previous article, and obtain 



U 



arbitrary^ 



= u, 



arbitrary • 



(8) 
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In Fig. 12 we show examples of unit cell control, where appropriately chosen 
parameters C^."^ and variations 5Ck achieve unitary transformations on the unit 
cell: the Toffoli-gate transformation (see Appendix two-qubit permutations 
Pij\a)i\b)j — \h)i\a)j (a, 6 = 0, 1), and the conditional phase shift employed in the 
quantum discrete Fourier transform (discussed in Sec. The transformation 
is achieved either directly (n* = 1) or by 8 repetitions in* — 8). The operators 
iJo, -Pw and Ps are chosen with arbitrary realistic values. We take Dyi = 
I?23 = OMQEu, Di3 = OME^, and T = 250h/Eu, where E^ ~ lO"!^ erg is the 
typical energy scale. For odd k we switch off the external electromagnetic field, 
^i;2;3 = 0, and tune the atomic excitation energies by the Stark field Es such 
that Ai;2;3 = (0.1; 0.11; 0.312)£:„. For even k we set Es = 0, that is Ai;2;3 = 0, 
and take Vi.2;3 = (0.3; 0.33; 0.24)Eu. 
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Figure 2: (a) Control parameters C° for the identity transformation /. Vari- 
ations SCk achieving the transformation U^/g, on the cell, with J/^ = ([/^/g)* 
equal to: (b) the permutation Upi2; (c) the permutation Up23; (d) the Toffoli- 
gate transformation Uxoff- (e) Variations SCk achieving the conditional phase 
shift 3(4)) = e'iq){—i(f)HB), at (j) — 7r/32, employed in the quantum discrete 
Fourier transform. 
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3 Completely Controlled Quantum Devices 



Once completely controlled unit cells have been constructed, a compound device 
can be assembled from such elements in an architecture which depends on the 
specific function it has to perform. Fig. |3| shows two possible arrangements 
of unit cells designed for two different purposes : the first one suites more the 
purpose of quantum computing, while the second is more useful for simulating 
lattice quantum field dynamics. 

The first device (Fig. EJa)) is organized in a tree-like structure. In this 
arrangement, the quantum state of one atom in each cell can be exchanged with 
the state of an atom in the closest parent joint of the tree. Even though the 
simplest way to make the exchange is to straightforwardly displace the atom 
to the parent joint, the exchange or transport of the state without moving the 
atom can be more practical. The tree-like architecture and the possibility to 
perform all the unitary transformations (including all the permutations) in each 
unit cell allow one to put together the states of any three two-level atoms of 
the device and make them interfere after at most s — 6 logg n state exchanges, 
by moving them toward the root of the tree to a common cell. Placing the new 
states back (if needed) requires the same number of inverse exchanges. This is 
a very modest number, s ^ 40, even for a rather large device of n ~ 10"^ with 
Hilbert space of iV = 2" ^ 10'^''° dimensions. Hence, all basic operations of 
quantum computation can be performed on any physical system composed of 
non-holonomic triads of two-level subsystems arranged in a tree-like structure, 
and each operation can be completed within 64xl6xl2x log3 n control intervals 
T. Note that the identity transformation should be applied to all other cells to 
preserve their states during the operation. 

The second arrangement of cells (Fig. |3Jb)) is mainly designed to emulate 
the dynamics of quantum fields on lattices. Of course, it can also perform 
general operations on any triad, but for a higher cost of s = 0(?i^/^). In this 
arrangement, after each control period of 64r the closest neighboring atoms 
are differently regrouped in triads (cells), with the original grouping repeating 
itself after three consecutive periods. Therefore, at each moment the change of 
the cell state depends on the states of the neighboring cells, as it should be in 
order to emulate the dynamics of the fields. Immediate analogy to the Ising 
model emerges when we restrict ourselves to small values of e for which terms of 
order are negligible, and then each T — 64r period plays the role of the time 
increment At = e. The evolution of such a device is determined by three sums 
of effective cell Hamiltonians, H^^^j: — J2q'^q,p^ O"^^ fo'" ^ach period p = 1,2,3, 
where Ti-q^p is the effective Hamiltonian of the qth cell at the pth period. 

We can cast the cell Hamiltonians to sums of tensor products of Pauli matri- 
ces CT^, where the Greek index a = x,y, z denotes the matrix type and the Latin 
index i specifies the two-level atom on which it acts. Since the cells are under 
complete control, the coefficients of this development can be made an arbitrary 
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Figure 3: Two possible arrangements of cells for special purpose devices: (a) 
tree-like structure for quantum computation; (b) planar lattice for simulating 
dynamics of quantum fields. The circled numbers denote the rank of joints of 
the tree (a) or specify the order in which atoms are grouped into triads (b). The 
arrows show state exchange to parent joints. 
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function of the time r, and hence the effective Hamiltonian reads 



(9) 



with impUcit summation over repeated indices, where («, j) and (ijj, A;) indicate 
pairs and triads of distinct atoms that are periodically grouped in a common cell. 
This Hamiltonian results in the evolution equation for the Heisenberg operators 



where the coefficients B,C are determined by A, B, C and the commutation 
relations of the Pauli matrices. By a proper choice of the coefficients A, B, C 
through the appropriate control sequences, one can simulate different linear and 
non-linear lattice models of quantum fields with time dependent parameters. 

4 Toy Device 

To conclude this paper, we now describe a toy device that can perform quan- 
tum computations on 9 qubits. An ensemble of 9 different Rydberg atoms, i. e. 
atoms of different elements or identical atoms which are excited to distinct pairs 
of Rydberg states, is placed in a magneto-optical trap at low temperature, as 
illustrated in Fig. 0] The best candidates for such a device are the long-living 
states corresponding to large angular momentum. By placing all the atoms in 
a static electric field, one lifts the degeneracy in the magnetic quantum number 
and performs tuning if needed. All the atoms experience the dipole-dipole inter- 
action Dij = didj{R~^) , where the cube of the inverse distance between atoms 
is averaged over their translational quantum states. Note, however, that only 
for almost resonant atoms this interaction is important. By a proper choice of 
the atomic states and the static field Es, we obtain three triads, p = 1,2, 3, each 
of which is composed of three almost resonant two-level atoms with transition 
frequencies centered on distinct frequencies ujp. For each triad p, the interac- 
tions Dij give the principal Hamiltonian, while a microwave field E^jj, at the 
frequency Up serves as a control perturbation. Transportation of the state of 
one atom in each triad to the parent joint can be performed by dipole or Raman 
TT transitions from the initial pair of Rydberg levels to a higher pair. With these 
higher pairs assumed nearly resonant with a frequency UJ4, atoms 3, 6 and 7 
form a higher-level triad — the parent joint of the first three triads — which is 
controlled by a fourth microwave field E^^ of frequency UJ4. 

As an application of our technique, we show how to perform the discrete 
Fourier transform modulo = 2^ = 512 on the toy device we have just pre- 
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Figure 4: A toy device, composed of 9 Rydberg atoms, which can perform 
quantum computations on 9 qubits. Each atom is a two-level system shown 
schematically by double orbits. Atoms of different triads arc excited to distinct 
pairs of Rydberg states. Each triad p is controlled by an external field of distinct 
frequency ujp. One atom in each triad can be excited to a pair of higher Rydberg 
states, thus forming a higher- level triad: (3,6,7). These excitations (depicted 
by arrows) correspond to state transportations. 
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sented. This operation is achieved by the following 9-qubit unitary transforma- 
tion 



JV-l 



Fn\x) = 



exp{2TTixy/N)\y), 



(11) 



v=o 



where \x) and \y) belong to the system computational basis, the states of which 
are defined by 

\x) = \xs)9 ■ ■ ■ |a;i)2|a;o)i, (12) 

with X = J2r=o ^r'^^ = 0, 1, . . . , iV — 1 {xr = 0, 1), and where | )i denotes the 
state of the ith atom — the ith qubit. The algorithm we employ to perform the 
Fourier transform is based on constructing the exponent in Eq. as 



exp(27ria;2//2^) = J| J]^ exp(i7r42/s/2'^ 



(13) 



r=0 s=0 



where xj, = a;8-r- We begin by reversing the order in which the bits of the input 
X are stored in our 9-qubit register, that is, we achieve the unitary transforma- 
tion 

|a;8)9 • ■ • |a;i)2|a;o)i ^ |a;o)9 . . . |x7)2|x8)i (14) 

by applying a sequence of state exchanges. Then we complete the transform 
in the following 9 steps: (i) We "split" the first qubit (the state of atom 1) by 
applying the unitary transformation 



A: 



1 



V2 ( 1 



1 1 



exp 



— ITT 

7f 



1 - \/2 1 

1 -1 - V2 



(15) 



which maps |0) ^(|0) + |1)) and |1) ^(|0) - |1)). Note that this would 
already complete the Fourier transform if we had only one qubit. (ii) Next, 
we apply to the first and second qubits the conditional phase shift |a)2|fe)i 



piTrah/2 



a)2|&)i (a, 6 = 0, 1), explicitly given by 



B21 — 



/ 1 

10 

1 

\ 



B(7r/2), 



(16) 



where B{(j)) is the unitary transformation 
i?((/)) = exp 



/ 





V 



\ 





-1 / 



(17) 



Then we "split" the second qubit by applying the transformation A. This 
accounts for the contribution of the second most significant bit of the input 
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X. (iii) Similarly, in steps i — 3, 4, . . . , 9 we apply the conditional phase shift 
\a)^\b)j ^ e*™''/2""'|a)i|6)j (a, 6 = 0,1), that is. 



/ 1 \ 

10 

1 

V e*'^/^'-^ J 



B[ir/T-^), (18) 



to each pair of qubits (z, j), j = 1, 2, .., i — 1, and then apply the transformation 
Ai = A to the ith qubit. Note that after the ith step the first i qubits store the 
Fourier transform of the i most significant bits of x. Hence, after the 9th step 
the Fourier transform is completed: 

= (IgSgg • • • Sgi) • • • {MB:i2B3i){A2B2l){Ai) . (19) 

Performing these operations implies also application of state exchanges when- 
ever one needs to transfer the states of atoms i and j to a common unit cell 
for processing. A list of control commands {5Ck sequences) corresponding to 
Eqs. H14() and H19|l can be written straightforwardly. 

5 Conclusion 

Non-holonomic control is very relevant in the context of quantum computation, 
since it allows one to perform any unitary evolution of the system chosen as 
a computer, or, in other terms, it is a way to implement any quantum gate. 
Moreover, direct control over the entire system is not always necessary in order 
to perform computations on the information it contains : indeed, a completely 
controlled quantum device can be constructed as the smart arrangement of 
universally controlled unit cells. Two examples of such architectures have been 
provided in this paper, which are built in accordance with the computational 
task they are supposed to carry out. Moreover, as an example, a toy device has 
been considered, which is able to perform a concrete computation on 9 qubits. 

A Toffoli gate 

The Toffoli-gate transformation is the unitary transformation on three qubits, 
UToff\x2)\xi)\xa) = |x2)|xi)|xoXOR(2;i ANDa;2)), (20) 
which corresponds to the three-bit classical logic gate, 

X2 X2 = X2 

xi x'l = Xi 

xo a;;, = a:oXOR(a;i ANDa;2), (21) 
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introduced by TofFoli as a universal gate for classical reversible computation. 
It acts as a permutation of the computational basis states, \x 
X = J2r=o 2;r2'" = 0, 1, . . . , 7, given by the unitary matrix 



/ 



^Toff 



This matrix can be presented as 

Uroff = expi-inHToff), 
with the (idempotent) Hcrmitian matrix 



Ti-Toff 



\ 










-1 
-1 1 



F2;fi)|xo) 



(22) 



(23) 



(24) 



In our control scheme the Toffoli-gate transformation can be effected on the unit 
cell by repeating 8 times the transformation U^/s = exp(— zttTYto/z/S) (e = tt), 
which is directly attainable: U{t = 64T) = C/^/g (see Fig.l^Id)). 
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